 plot.pacific.skj <-
function(plot_title="") {
eez <- read.table("L:\\alb\\2008\\Pago\\eznew2.txt")
#source("R addpie.r")
plot(1,1, yaxt="n", xaxt="n", type="n", xlim=c(110,220), ylim=c(-40,45), ylab="", xlab="", bg="lightblue")
polygon(c(100, 230, 230, 100), c(-45, -45, 50, 50), col="lightblue")
polygon(eez[,1], eez[,2], lwd=1, col="white")
lines(eez[,1], eez[,2], lwd=1, col="slate grey")
map('world2Hires',  yaxt="n", xaxt="n", add=T, resolution=1)
map('world2Hires',  region = c("USA","Hawaii","Mexico","Japan","China","South Korea","North Korea","Philippines","Vietnam","Laos","Taiwan","Fiji", "Vanuatu", "Malaysia", "Australia", "New Zealand", "Indonesia", "New Caledonia", "Papua New Guinea", "Solomon Islands"), fill=T, add=T, yaxt="n", xaxt="n", col="black", density=50)
box(lwd=3)
lines(c(140,140,210), c(25,45,45), lwd=2,lty=1)
lines(c(120,120,140), c(25,35,35), lwd=2,lty=1)
lines(c(165,165), c(25,45), lwd=2,lty=1)
lines(c(115,210), c(25,25), lwd=2,lty=1)
lines(c(115,210), c(15,15), lwd=2,lty=1)
lines(c(115,115,142.5,142.5,210), c(25,-10,-10,-20,-20), lwd=2,lty=1)
lines(c(165,165), c(15,-20), lwd=2,lty=1)
lines(c(210,210), c(45,-20), lwd=2,lty=1)
text(138,27,labels=1,font=2,cex=1.2)
text(138,27,labels=1,font=2,cex=1.2)
text(163,27,labels=2,font=2,cex=1.2)
text(208,27,labels=3,font=2,cex=1.2)
text(208,18,labels=4,font=2,cex=1.2)
text(163,-18,labels=5,font=2,cex=1.2)
text(208,-18,labels=6,font=2,cex=1.2)
axis(1, at=c(100,110,120,130,140,150,160,170,180,190,200,210,220,230,240,250), labels=F)
axis(2, at=c(-40,-35,-30,-25,-20,-15,-10,-5,0,5,10,15,20,25,30,35,40), labels=F)
axis(2, at=c(-35,-25,-15, -5,5,15,25,35), labels=c("35S","25S", "15S","5S","5N","15N","25N","35N"), cex.axis=0.75)
axis(1, at=c(120,140,160,180,200,220,240), labels=c("120E", "140E", "160E", "180", "160W", "140W", "120W"), cex.axis=0.75)
mtext(side=3, line=0.5, plot_title)
}
